Percolation in three-dimensional random field Ising magnets 
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The structure of the three-dimensional random field Ising magnet is studied by ground state 
calculations. We investigate the percolation of the minority spin orientation in the paramagnetic 
phase above the bulk phase transition, located at [A/J] c ~ 2.27, where A is the standard deviation 
of the Gaussian random fields (J = 1). With an external field H there is a disorder strength 
dependent critical field ±_ff c (A) for the down (or up) spin spanning. The percolation transition 



where A,- 



2.43 ± 0.01 and 



PACS 



is in the standard percolation universality class. H c ~ (A — 71 
S = 1.31 ± 0.03, implying a critical line for A c < A < A p . When, with zero external field, A is 
decreased from a large value there is a transition from the simultaneous up and down spin spanning, 
with probability II n = 1.00 to II U = 0. This is located at A = 2.32 ± 0.01, i.e., above A c . The 
spanning cluster has the fractal dimension of standard percolation Df — 2.53 at H — H C (A). We 
provide evidence that this is asymptotically true even at H = for A c < A < A p beyond a crossover 
scale that diverges as A c is approached from above. Percolation implies extra finite size effects in 
the ground states of the 3D RFIM. 
0,5.50.+q, 75.60. Ch, 75.50.Lk, 64.6O.Ak 



I. INTRODUCTION 



The random-field Ising model (RFIM) is one of the 
most basic models for random systems fl|-|3J . Its beauty 
is that the mixture of random fields and the standard 
Ising model creates rich physics and leaves many still 
un-answered problems. By now it is known that three 
dimensions (3D) is the corner stone of the model, since it 
presents a phase transition where the randomness proves 
to be a relevant perturbation to the pure 3D Ising model. 
For the last fifteen years, since the method-wise seminal 
paper by Ogielski Q, the studies of the transition have 
centered around zero-temperature ground state computa- 
tions because the temperature is due to renormalization 
group arguments believed to be an (perhaps dangerously 
so) irrelevant variable. 

Many such works exist so far, the most recent and 
comprehensive being due to Middleton and Fisher Q . In 
spite of all the effort many uncertainties remain concern- 
ing the nature of the phase transition. The question is if 
the transition is of the second-order, of traditional first- 
order type, or finally some other kind of discontinuous 
transition. The order-parameter exponent, 3, may have 
a finite value or it can be equal to zero lq--|l2]|. Its very 
small value makes it unlikely that insight is obtained in 
the near future, in spite of the fact that the optimization 
algorithms used can at best scale almost linearly with the 
number of spins in the system. Moreover a controversy 
exists with regards to the role of disorder: the available 
simulations are not able to settle the question whether 
the critical exponents depend on the particular choice of 
the distribution for the random fields, analogously to the 
mean-field theory of the RFIM where binary (±/i) disor- 
der results in a first-order transition and Gaussian (see 




below) in a second-order one 113 ]. 

In this paper we focus on a novel aspect of the three- 
dimensional RFIM, namely percolation jl4j . The goal is 
to explore percolation critical phenomcnona in the 3D 
RFIM. The work is an extension to our studies of perco- 
lation in two-dimensional RFIM In the traditional 
3D Ising model, without disorder, the percolation behav- 
ior in an applied field and its consequences, as whether 
the phase transition critical exponents would be affected 
by the percolation criticality, are known since long as the 
"Trouble with Kertesz" This problem was solved 

by Wan, 
Klein 



by studying Fortuin-Kasteleyn or Coniglio- 
clusters using so called ghost-spins. In the 
RFIM the situation is different in that at small temper- 
atures one has a non-zero spin-spin overlap q with the 
ground state: thus the existence of a ground state perco- 
lation transition (even without an external field) implies 
measurable consequences even at finite temperatures. It 
also complicates the phase diagram by its existence. 

There is one fundamental difference between two and 
three dimensions (besides the fact that there is no phase 
transition in two dimensions, and hence there systems 
are always paramagnetic). In two dimensional square lat- 
tices the critical percolation site-occupation probability 
is 0.592746, i.e., above one half, and in three dimensional 
cubic lattices well below one half, 0.3116. Therefore in 
three dimensions, deep in the paramagnetic phase, both 
of the spin orientations should span the system (this has 
been noted by Esser et al. to be true for the RFIM, see 
pf| ) . Thus introducing an external field in paramagnetic 
systems leads in two dimensions to the percolation of the 
spin direction parallel with the external field. In three di- 
mensions, on the other hand, the external field destroys 
the spanning property of the spin orientation opposite to 
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the external field. 

Consequences of the percolation type of order at the 
paramagnetic phase are many-fold. There are experi- 
mentally accessible random field magnets, so called di- 
luted antiferromagnets in an external field (DAFF) pl| ] 
in which the percolation order could be seen, should it 
exist for zero external fields. It is already known that 
the percolation of the diluted atoms has a strong contri- 
bution to the behavior of the structure factor line-shapes 
of the 3D DAFF ||-|||. Near the thermodynamical 
phase transition point the universality class of the tran- 
sition is determined by several exponents, among them by 
the correlation length exponent, (if the transition is con- 
tinuous). The critical percolation phenomenon near the 
thermodynamical phase transition point may contribute 
there and introduce extra corrections, which have to be 
taken into account when the thermodynamical correla- 
tion length exponent is determined. 

This paper is organized so that it starts with an in- 
troduction of the random field Ising model in the next 
section. Also the numerical method solving exactly the 
ground states is introduced. In Section [II the perco- 
lation phenomenon is studied, with a non-zero external 
field. The universality class of the percolation behav- 
ior is determined and the dependence of the critical ex- 
ternal field on the random field strength is investigated. 
Section [Fv| concentrates on the percolation phenomenon 
without an external field and compares it with the cases 
when the external field is applied. The pro perties of the 
spanning cluster are studied in Section [Vj Implications 
of the percolation to the phase diagram are discussed, 
together with the conclusions in Section |v|. 



II. RANDOM FIELD ISING MODEL AND 
NUMERICAL METHOD 

The random field Ising model is defined by its energy 
Hamiltonian 



(v) 1 



(1) 



where J > (throughout this paper we set J = 1, 
since the relevant value is its ratio with the random 
field strength) is the coupling constant between nearest- 
neighbor spins Si and Sj. We use here cubic lattices. H 
is a constant external field, which if non-zero is assigned 
to all of the spins, and hi is the random field, acting on 
each spin Si- We concentrate only on a Gaussian distri- 
bution for the random field values 



P(fH) = 



1 



2ttA 



exp 



2 I A 



(2) 



with the disorder strength given by A (in this paper A 
actually denotes the ratio between disorder strength and 



the coupling constant), the standard deviation of the dis- 
tribution. The arguments presented in this paper could 
be extended to other lattices and other distributions, e.g. 
uniform and bimodal, too. However, discrete distribu- 
tions, such as the bimodal one, suffer from degeneracies, 
and when calculating thermodynamical quantities extra 
averaging, over the degeneracies, has to be done when 
using discrete distributions |p5| , p6| . 

To find the ground state structure of the RFIM means 
that the Hamiltonian (|l|) is minimized, in which case the 
positive ferromagnetic coupling constants prefer to have 
all the spins aligned to the same direction. On the other 
hand the random field contribution is to have the spins to 
be parallel with the local field, and thus has a paramag- 
netic effect. This competition of ferromagnetic and para- 
magnetic effects leads to a complicated energy landscape 
and the finding the ground state becomes a global opti- 
mization problem. An interesting detail of the RFIM is 
that for H — it has an experimental realization as a di- 
luted antiferromagnet in a field. By gauge-transforming 
the Hamiltonian of DAFF 



H 



(ij) 



(3) 



where the coupling constants J < 0, e$ is the occupation 
probability of a spin Si , and B is now a uniform external 
field, one gets the Hamiltonian of RFIM (|l|) with H = 
p7| , p8 21 1 . The ferromagnetic order in the RFIM corre- 
sponds to antiferromagnetic order in the DAFF, natu- 
rally. 

For the numerical calculations a graph-theoretical 
combinatorial optimization algorithm developed in com- 
puter science has been used. The Hamiltonian ([!]) is 
transformed to a random flow graph widely used in com- 
puter science with two extra sites: the source and the 
sink. The positive field values hi correspond flow capaci- 
ties Cn connected to the sink (t) from a spin Si, similarly 
the negative fields with a s are connected to the source 
(s), and the coupling constants 2J^ 
spins correspond flow capacities = 
to its neighboring one Sj In the case the exter- 

nal field is applied, only the local sum of fields, H + hi, 
is added to a spin toward the direction it is positive. 
The algorithms, namely maximum-flow minimum-cut al- 
gorithms, enable us to find the bottleneck, which restricts 
the amount of the flow which is possible to get from the 
source to the sink through the capacities, of such a ran- 
dom graph. This bottleneck, path P which divides the 
system in two parts: sites connected to the sink and sites 
connected to the source, is the global minimum cut of the 
graph and the sum of the capacities belonging to the cut 
Y]p Cij equals the maximum flow, and is smaller than of 
any other path cutting the system. The value of the max- 
imum flow gives the total minimum energy of the system 
and the minimum cut defines the ground state structure 
of the system, so that all the spins in the source side of 



= between the 
Cji from a site Si 
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the cut are the spins pointing down, and the spins in the 
sink side of the cut point up. The maximum flow algo- 
rithms can be proven to give the exact minimum cut of 
all the random graphs, in which the capacities are pos- 
itive and with a single source and sink |30). We have 
used a sophisticated method for solving the maximum 
flow - minimum cut problem called push-and-relabel by 
Goldberg and Tarjan |3l| , which we have optimized for 
our purposes. It scales almost linearly, C(n 12 ), with 
the number of spins and gives the ground state in about 
minute for a million of spins in a workstation. 

We have used periodic boundary conditions in all of 
the cases. Also the percolation is tested in the periodical 
or cylindrical way, i.e., a cluster has to meet itself when 
crossing a boundary in order to span a system. Find- 
ing the spanning cluster has been done using the usual 
Hoshen-Kopelman algorithm |3^] . 



III. PERCOLATION WITH AN EXTERNAL 
FIELD 

As a start of the percolation studies of the 3D RFIM we 
draw in Fig. |l|(a) the spanning probabilities of down spins 
II| with respect to the uniform external field H pointing 
up for several system sizes L and for a fixed random field 
strength A = 3.5. The curves look rather similar to stan- 
dard percolation, except that in site percolation the sys- 
tems span at high occupation probability limit, and here 
the down spins do not span, when the external positive 
field has a large value, and thus the step in the spanning 
probability is inverse compared to the one in the occupa- 
tion percolation. It is interesting to note also, that since 
we are using periodic boundary conditions in all of the 
directions, also for spanning, the n^(L) -lines for various 
system sizes cross at rather low LT^ values. This is the 
case for the other A, too. Similar boundary condition 
dependent behavior have been seen in the standard per- 
colation, too ]3^-|35|]. When we take the crossing points 
H C (L) of the spanning probability curves with fixed span- 
ning probability values Hi = 0.4, 0.5, 0.6, 0.7, and 0.8, 
for each systems size L, we get an estimate for the critical 
external field H c using finite size scaling, see Fig. |l|(b). 
There we have attempted with success to find the value 
for H c using the standard short-range correlated 3D per- 
colation correlation length exponent v — 0.88 |L4|]. Using 
the estimated H c = 0.461 ± 0.001 for A = 3.5 we show a 
data-collapse of IT^ versus (H — H c )/L~ 1 / 1 ' in Fig. |l|(c), 
which confirms the estimates of H c and v — 0.88. We 
get similar data-collapses for various other random field 
strength values A as well. 

Considering the percolation and critical external field 
with respect to the random field strength, there is an ob- 
vious constraint in the phase diagram H vs. A. Below 
the phase transition critical point, A c ~ 2.27 |§,0J|], only 
one of the spin orientations may span a system, since in 
a ferromagnetic system the magnetization has a finite, 



positive or negative, value and thus there can not be a 
massive percolation cluster of the opposite spin direc- 
tion. Since the earlier studies of the phase transition at 
3D RFIM ||-||,(i]]jl2) have shown that the order param- 
eter exponent (3 has a value close to zero, if not zero, the 
transition is sharp and therefore the simultaneous perco- 
lation of the both (up and down) spin directions should 
vanish or have vanished at A c when approaching from 
above. The question now remains, whether this takes 
place exactly at the phase transition point, so that the 
critical points would coincide, or for a A p > A c . In the 
latter case it is also of interest what happens for H = 
between the critical points, on the line A c < A < A p . 
We now propose a phase diagram, Fig. Q, for the perco- 
lation phenomenon, and ask at which value the dashed 
lines in the diagram meet. Above we showed that in 
the direction of the vertical arrow at H > the universal 
standard percolation correlation length exponent is valid. 
What about at the vertical arrow, what are the critical 
exponents there? 

To answer the question how the percolation critical ex- 
ternal field H c behaves with respect to the random field 
strength A, we have attempted a critical type of scaling 
using the calculated H C (A) for various A =2.5, 2.6, 2.75, 
3.0, 3.25, 3.5, 4.0, and 4.5. We have been able to use the 
Ansatz 



H c ~ (A - A p ) s 



(4) 



where S = 1.31 ± 0.03 by assuming A p = 2.43, see 
Fig. ||(a). In Fig. ||(b) on the other hand we have plot- 
ted the calculated A values versus the scaled critical ex- 
ternal field [-ffc(A)] 1 / 1 - 31 and it gives the estimate for 
A p = 2.43 ± 0.01. This indicates that the percolation 
probability lines for up and down spins to lose their span- 
ning property meet at A p = 2.43 ± 0.01. Note, that 
our studies in two-dimensional RFIM gave the values 
A p = 1.65 ± 0.05 and S = 2.05 ± 0.10 for systems to 
span, not to lose the spanning property as here [ fL5"| . We 
also tested various exponential scaling assumptions for 
the H C (A) scaling, but none of them worked. However, 
here we know, that H c has to vanish at some finite A p 
value, which is greater than or equal to A c . 

We have also calculated the order parameter of the 
percolation, the probability that a down spin belongs to 
the down-spin spanning cluster P^. Using the scaling 
for the correlation length 

£,perc ~ \H — H c \ , (5) 

and for the order parameter, when L < £, pe rc, 

P 00 (H)r,(H c -Hf, (6) 



we get the limiting behaviors, 
Poo{H,L)~{f_%J v H ^ 



^* £,perci 



(7) 
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and thus the scaling behavior for the order parameter 
becomes 

P OB {H t L)~L-V''F \ {Hc ~ L HrV 

Note, that here and later in this article (3 denotes the per- 
colation order parameter exponent as opposed to the bulk 
phase transition order parameter exponent discussed ear- 
lier in this paper. We have done successful data-collapses, 
i.e., plotted the scaling function /, for various A using 
the standard 3D short-range correlated percolation expo- 
nents /? = 0.41 and v — 0.88, of which the case A = 4.5 
with H c = 1.0441 is shown in Fig. ||. Note, that only the 
left part (below zero) of the scaling function is shown, 
since /^(if, L) is limited between [0,1]. When one di- 
vides it by L~PI V the part where non-scaled Poc(H : L) 
had a value of unity the scaled P 00 (7f, L)/L~P/ V satu- 
rates at different values depending on L. One can eas- 
ily see, that the smallest system size L 3 = 8 3 does not 
scale (the rest are scattered around each other and do 
not have any trend). We believe that this is due to an 
intrinsic length scale over which the spins are correlated, 
and which depends on the random field strength value. 
This will be discussed in more detail in Section when 
the scaling of the spanning cluster is studied. 

Hence, we conclude that the percolation transition for 
a fixed A versus the external field H is in the standard 3D 
short-range correlated percolation universality class |L4j . 
This is confirmed by the fractal dimension of the span- 
ning cluster, too, as discussed below. The fact that the 
critical behavior of the percolation with respect to the 
external field belongs to the standard short-range corre- 
lated percolation universality class is not surprising, since 
the strong disorder limit can be seen to be related with 
the site percolation problem and that e.g. the positive 
external field decreases the number of the occupied down 
spins. Also other exponents could be measured, as 7 for 
the average size (s) of the clusters, and a and r for the 
cluster size distribution as well as the fractal dimension 
of the backbone of the spanning cluster, the fractal di- 
mension of the chemical distance, the hull exponent etc. 

IV. PERCOLATION AT H = 

In the previous section we learned that the dashed 
lines at the phase diagram, Fig. |^, meet at the value 
A p = 2.43 ± 0.01, which is well above the phase transi- 
tion critical point A c = 2.27. This raises the question, 
how this is seen, when the external field H = and what 
happens between A c and A p . Thus we study the phase 
diagram in the direction of the horizontal arrow in Fig. |^. 
There are two strategies for this that we employ sepa- 
rately to evaluate their advantages and disadvantages. 



That is, one can take the A p to be a priori the same for 
all Il-fi, the probability for simultaneous spanning of up 
and down spins. Or then this can be let to vary with 
IIf|, as in two-dimensions p5| . 

In Fig. ||(a) we have plotted the probability for simul- 
taneous spanning of up and down spins II ^ as a function 
of the Gaussian random field strength A for various sys- 
tem sizes L 3 = 8 3 , 15 3 , 30 3 , 50 3 , 90 3 , and 120 3 . This 
case now resembles the standard occupation percolation 
in the sense, that the step in the percolation probability is 
from a low value to a large value when A is increased. By 
estimating that the A p ^h=o at the thermodynamic limit 
has a value of 2.32 using fixed rifj, = 0.2, 0.4, 0.6, and 
0.8 for the A Pi h=o(L) we find that the effective v gets a 
value of 0.97 ± 0.05 when approaching the critical point 
in this direction, see Fig. |B|(b) . On the other hand assum- 
ing that the v = 1.0 the A p ,h=o becomes 2.32 ±0.01, see 
Fig. |B|(c) . These plots show that the estimates should be 
correct. However, the data-collapse, Fig. |](d), using the 
estimates above could be better. Obviously the smallest 
system size, L 3 — 8 3 , does not scale. 

There are a couple of points one should note from the 
scaling. Firstly the estimate for the Ap i ^f = o = 2.32 ±0.01 
is still above the phase transition point A c = 2.27. An- 
other point is that Ap ,h=o is reasonably far away from 
Ap = 2.43 ± 0.01 [note, that the error-bar in the fi- 
nite field case is the error-bar of the least-squares fit in 
Fig. ||(b) and does not take into account other sources for 
the error, e.g., the error of <5, statistics etc., and thus is a 
lower limit]. The third point is that H-\i = 0.0 at A„h=o 
and n x = 0.25 at H C (A) [for A = 3.5 see Fig. 0(c)]. 
Our take on the two different estimates is that they are 
compatible with the following scenario. For A values 
that are slightly below 2.43 one can have only one criti- 
cal spanning cluster, and the probability for this is then 
n^, about 0.25. The both orientations do span simulta- 
neously, as they can do for all A- values above A c , but 
they should not be both critical, unless one decreases the 
disorder strength further. 

For the estimate of the correlation length exponent, 
deviations from normal percolation are seen since v = 
1.00±0.05 instead of v = 0.88. In our opinion this reflects 
the fact that for H 7^ the correlations from the proxim- 
ity of A c are negligible, whereas here the spin-spin corre- 
lations change with system size. The correlation length 
exponent is higher than that for percolation, so clear-cut 
percolation scaling can not be expected. Differences be- 
tween the H = and H =^ -cases were found also in the 
two dimensional case Jig] . Note, that in two dimension, 
the exponent was dependent on the spanning probability 
and the standard correlation length exponent was found 
where the spanning probability for either of the spin di- 
rections to span /| had a non-zero value (remember, 
that in two dimensional square lattices without an exter- 
nal field at large A neither of the spin directions span, 
and with small A either of them start to span). 
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Here we tried, as in two-dimensions, to do fits using 
several criteria for IT n = [0.05, 0.15, 0.20, . . . , 0.95] and 
letting both A p and v vary depending on H-|-j_. Indeed, 
we obtained monotonous behaviors depending on IL^ for 
both v and A p . However, this may just reflect how finite 
size effects depend on the criterion. It is anyhow worth of 
noting that for Il-fi approaching zero, A p gets also closer 
and closer to 2.27, i.e., the accepted value for the phase 
transition point A c . Moreover the correlation length ex- 
ponent moves towards v = 1.3±0.1, in the neighborhood 
of the phase transition correlation length exponents re- 
ported in the literature p|,p|,pl]| . Similarly if II n is let 
to approach unity, A p closes on the value A p = 2.43 
obtained above, in the finite field case. This behavior 
may be just coincidence, or related to the (A-dependent) 
correlations in the system, to how they change the uni- 
versality class of percolation in the vicinity of the phase 
transition. We return to this in the conclusions, in Sec- 



tion VI 



Hence we have shown that at large A both of the spin 
directions span simultaneously, and by decreasing ran- 
dom field strength we find a critical A p .# = o, which is 
above the phase transition point A c , and below which 
there is no simultaneous spanning. Therefore we con- 
clude, that in the whole regime A c < A < A p h=o there 
is geometrical criticality in 3D RF magnets, since always 
only either of the spin directions spans the system. How- 
ever, the spanning cluster cannot be massive there, i.e., 
scale with the Euclidean dimension (d = 3), the system 
still being paramagnetic, but has to be a fractal. The 
scaling of the spanning cluster is studied in the next sec- 
tion and the implications of the critical region in Sec- 



tion VI 



V. THE SPANNING CLUSTER 

In Fig. ||(a) we have plotted the mass of the spanning 
cluster of down spins with respect to the system size at 
He (A) > for four random field strength values A = 
2.75, 3.0, 3.5, and 4.0 up to the system size L 3 = 120 3 . 
As a guide to the eye the fractal dimension Df — 2.53 
of the standard percolation is drawn in the figure and 
the systems can be seen asymptotically approaching the 
same scaling. However, there are obvious finite size ef- 
fects, which depend on A. We have estimated roughly 
the crossover system sizes for the systems to reach the 
correct scaling, L x ~ 30, 20, 10, and 5 for A = 2.75, 3.0, 
3.5, and 4.0, respectively. This hints about an exponen- 
tial scaling with a slope of — 1.42 ± 0.03 for the crossover 
length scale, see closed diamonds in the inset of Fig. ||(a). 
The above scaling predicts for A = 4.5 L x ~ 3, smaller 
than L = 8 (in Fig. ||, this size does not scale) but note 
that the prefactors of the scaling behaviors need not to be 
the same. In Fig. |^(b) we have drawn for three A < A p , 
i.e., A =2.35, 2.38, and 2.45 (which is so close to A p 



that its H c is practically zero with respect to the numer- 
ical precision, 10~ 3 ) at H = the scaling of the mass of 
the spanning cluster of either of the spin orientations up 
to system size L 3 = 120 3 . There one can see that the 
fractal dimension Df = 2.53 of the standard percolation 
is asymptotically met, too, but at much larger system 
sizes. Here we have estimated the crossover system sizes 
L x ~ 80, 60, and 50, for A =2.35, 2.38, and 2.45, respec- 
tively. They are plotted as open circles in the inset of 
Fig. ra(a) and are obviously diverging from the exponen- 
tial behavior mentioned above when approaching phase 
transition A c . These large values for L x do not leave 
much room for the asymptotic scaling, since it is difficult 
to go above L 3 = 120 3 . However, the crossover is visible. 
There is one another thing one notices from Figs. |^(a) 
and (b). In the case we plot the mass of the spanning 
cluster of the down spins in A > A p and H C (A) > 
the crossover is from a smaller slope to the asymptotic 
Df = 2.53 one. In the case A < A p the crossover is from 
the Euclidian dimension (slope of three, i.e., effective fer- 
romagnetism) to the asymptotic Df — 2.53. There it is 
obviously affected by the vicinity of the phase transition 
point. 



VI. DISCUSSION AND CONCLUSIONS 



In this paper we have studied the character of the 
ground state of the three-dimensional random field Ising 
magnet in, mostly, the paramagnetic phase. A geometri- 
cal critical phenomenon exists in these systems: for cubic 
lattices in ordinary percolation both occupied and unoc- 
cupied sites span the systems, when the occupation prob- 
ability is one half. In the RFIM this corresponds to the 
case with a high random field strength value, without an 
external field. When an external field is applied and the 
random field strength decreased, a percolation transition, 
for the other spin orientation to lose the spanning prop- 
erty, can be seen. The transition is shown to be in the 
standard 3D short-range correlated percolation univer- 
sality class, when studied as a function of the external 
field. Hence, the correlations in the three-dimensional 
random field Ising magnets are only of finite extent as 
could be expected in this region of the bulk phase dia- 
gram. Based on our numerical results both the critical 
points ±H C (A) approach when A is decreased, and fi- 
nally meet at a A p ~ 2.43 > A c , at which H c = 0. 
When the percolation transition is studied without an ex- 
ternal field and tuning the random field strength similar 
behavior is found, i.e. signatures of a percolation line (a 
A p > A c ). This might cause puzzling consequences when 
studying the character of the ground states, because the 
percolation correlations may influence the magnetization 
correlation length. 

The major theoretical implications have to do with the 
phase transition. Note that earlier groundstate studies of 
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the domain structure implied that there is only a "one- 
domain state" below the critical field, and a "two-domain 
state" in the paramagnetic phase (extending down from 
high disorder values) fpo[ . If the transition is first-order, 
then one expects the percolation properties of the param- 
agnetic phase to be discontinuous in the thermodynamic 
limit. If the transition is second-order, then one may 
ask what is the correct way to link the presence of the 
percolation transition to the critical phase? At A c , one 
expects that the spin-spin correlations show power-law 
correlations. For a normal percolation transition, these 
are (as in the disordered phase in general) of short-range 
character. There is a divergent length scale as the tran- 
sition is approached from the paramagnetic phase, below 
which the spin-spin correlations matter and the scaling 
of the spanning clusters is volume-like. 

Assume that the properties of the largest cluster are 
governed by the power-law correlations. An old result 
by Weinrib gives a Harris' criterion for this approach, to 
check how this would change its structure from ordinary 
percolation (3^] . If the site occupation probability corre- 
lations decay as r~ a , one has that the decay is relevant 
if avoid — 2 < —> v new = 2 /a, where now v i^ — 0.88 for 
3D site percolation. One gets a critical decay exponent 
a c = 2.27, much larger than that found by Middelton and 
Fisher j^], which is very close to zero. An application 
of the theory of correlated percolation would thus im- 
ply that the spin-spin correlations at A c are relevant for 
percolation. They would change the universality class, 
of percolation, in a way that would reflect such corre- 
lations. This conclusion should be taken with plenty of 
salt, obviously. 

One should note also that although this study was done 
using cubic lattices it can be extended to other lattices, 
too, since all the common three dimensional lattices have 
p c < 0.5. Thus the transition from the both spin orienta- 
tions spanning phase to only one spin orientation span- 
ning phase should exists. In the case of diluted antiferro- 
magnets the percolation is already seen as percolation of 
diluted spins. The implication of this paper is that the in- 
fluence of percolation is even more rich. Lately there has 
been interest in studying domain walls and excitations 
in RF magnets. In both cases the underlying percolation 
criticality should affect the structure of the clusters that 
result from varying the boundary conditions. 
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FIG. 1. (a) The spanning probabilities ol minority down 
spins II | as a function of upward external field H for A = 3.5 
with L 3 6 [8 3 — 90 s ] . The number of realizations varies be- 
tween 5000 realizations for L = 8 and 200 for L — 90. (b) 
The finite size scaling of the fields H C (L), which are from the 
crossing points of the spanning probability curves with the 
horizontal lines in (a), leading to the estimate of the critical 
H c = 0.461 ± 0.001 using L~ 1/l/ , v = 0.88. The error-bars in 
the labels of the figure for different H c are the errors of the 
least-squares fits, (c) The data-collapse of different system 
sizes with the corresponding critical H c = 0.4608. 
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FIG. 6. (a) The average mass of spanning cluster of down 
spins for random field strength values A = 2.75, 3.0, 3.5, and 
4.0 at the critical positive external field value H C (A) [see the 
values from Fig. [](b)]. The 3D percolation fractal dimen- 
sion Df = 2.53 is indicated with solid lines. In the inset a 
crossover length scale at which system size the asymptotic 
behavior is met for each random field strength is plotted as 
closed diamonds. The least-squares fit estimates an exponen- 
tial behavior with a slope of —1.42 ± 0.03. (b) The aver- 
age mass of spanning clusters of either spin orientations for 
random field strength values A = 2.35, 2.38, and 2.45 when 
H = 0. The solid lines are the least-squares fits to the data 
with the slopes indicated in the labels. The dotted line with 
a slope of Df — 2.53 and the dashed line with a slope of d = 3 
are guides to eye. The estimated crossover length scales are 
plotted in the inset of Fig. E3(a) as open circles. 
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